Plasmoid and Kelvin-Helmholtz instabilities in Sweet-Parker current sheets 
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A 2D linear theory ol the instability of Sweet-Parker (SP) current sheets is developed in the 
framework of Reduced MHD. A local analysis is performed taking into account the dependence of 
a generic equilibrium profile on the outflow coordinate. The plasmoid instability [Loureiro et al, 
Phys. Plasmas 14, 100703 (2007)] is recovered, i.e., current sheets are unstable to the formation of 
a large- wave- number chain of plasmoids (fc max ics ~ S 3//8 , where fc max is the wave- number of fastest 
growing mode, S — LcsVa/ti is the Lundquist number, Lcs is the length of the sheet, Va is the 
Alfven speed and t] is the plasma resistivity), which grows super- Alfvenically fast (7maxi~A ~ S 1 ^ 4 , 
where 7 max is the maximum growth rate, and ta = Lcs/Va)- For typical background profiles, 
the growth rate and the wave-number are found to increase in the outflow direction. This is due 
to the presence of another mode, the Kelvin-Helmholtz (KH) instability, which is triggered at 
the periphery of the layer, where the outflow velocity exceeds the Alfven speed associated with 
the upstream magnetic field. The KH instability grows even faster than the plasmoid instability, 
Lcs ~ S 1 / 2 . The effect of viscosity (u) on the plasmoid instability is also addressed. 
In the limit of large magnetic Prandtl numbers, Pm = u/n, it is found that 7 max ~ 5 ,1//4 Pm~ 5//8 and 
fcmaxics ~ S 3 ' 8 Pm~ 3 ' 16 , leading to the prediction that the critical Lundquist number for plasmoid 
instability in the Pm 3> 1 regime is S C vit ~ 10 4 Pm 1 ' 2 . These results are verified via direct numerical 
simulation of the linearized equations, using a new, analytical 2D SP equilibrium solution. 

PACS numbers: 52.35.Vd, 52.35.Py, 94.30.cp, 96.60.lv 



I. INTRODUCTION 

Magnetic reconnection (THl] is a ubiquitous plasma 
physics phenomenon, characterized by the rapid recon- 
figuration of the magnetic-field topology. Solar flares Q 
and magnetospheric substorms @ are two prominent ex- 
amples of events where reconnection plays a key role. 
Plasma dynamics in many laboratory experiments is also 
critically determined by magnetic reconnection; exam- 
ples are the sawtooth [tJ and the tearing instabilities 
in magnetic-confincment-fusion devices, or the reconnec- 
tion of high- energy- density, laser-produced plasma bub- 
bles [IM|. 

Along with fast reconnection rates, many observa- 
tions [13j of magnetic reconnection phenomena display 
one intriguing feature: the formation, and subsequent 
ejection from the current sheet, of coherent secondary 
structures, often referred to as plasmoids (also known 
as blobs, flux ropes or secondary magnetic islands). 
There is abundant direct evidence for the presence of 
these structures in the Earth's magnetotail (14H16| and 
in solar flares p7l - l22j . In magnetic confinement fu- 
sion devices, plasmoid generation seems to be less cer- 
tain, though there are reports of the observation of sec- 
ondary magnetic structures correlated to m/n = 1/1 and 
m/n = 2/1 magnetic islands on the TEXTOR [$1 [H[ 
and JET [H[ tokamaks. On TEXTOR, high-resolution 
measurements of electron temperature fluctuations show 
structures which hint at plasmoid formation during saw- 



tooth crashes j26H28| . Finally, recent laser-plasma exper- 
iments where reconnection is conjectured to occur also 
show evidence for plasmoid formation [ll|, [l2[ • 

Direct numerical simulations of reconnection processes 
concur with observations in displaying ubiquitous evi- 
dence for plasmoid formation. Plasmoids have been re- 
ported in numerical simulations using various physical 
models, ranging from kinetic |29|-[33l to Hall-MHD [34|, 
[HI and to single fluid MHD [30l - l49j ] . Plasmoid formation 
has also been reported in numerical simulations of recon- 
nection in relativistic plasmas, both resistive f50j and ki- 
netic (5ll - [54j . Numerical studies tailored to address spe- 
cific reconnection contexts such as the solar corona [55l - 
[HH, the Earth's magnetotail [59l [60j . magnetic young 
stellar objects [6ll - l63| . fusion experiments [36l l39j and 
laser- plasma interactions [64[ , though different from each 
other in a number of details, again all appear to agree on 
the basic fact that plasmoids are generated in reconnect- 
ing current sheets. 

The plasmoid dynamics inferred from observations 
and seen in numerical simulations strongly suggest the 
very opposite of the laminar, steady-state reconnec- 
tion scenarios that have dominated the field for much 
of its history [the Sweet-Parker (SP) |(35|, [6(| and the 
Petschek [631 models, and, more recently, the Hall re- 
connection paradigm (f38j]. Magnetic reconnection in 
the presence of plasmoids appears to be a highly time- 
dependent, bursty process, which can only be described 
in a statistical manner (49|, [HH, l(39l - [72i . Furthermore, in 
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addition to their key role in setting the reconnection rate 
in both laminar 0, HI, HE III \MM, H HI 111 and tur- 
bulent [HI, EH plasmas, there is numerical and observa- 
tional evidence that plasmoids may be critical in explain- 
ing electron acceleration in reconnection sites pjl I73l - l75j ] . 

In a previous paper [76| (henceforth referred to as Pa- 
per I) we attempted to understand the origin of plas- 
moid formation in reconnection sites by analysing the 
linear stability of large-aspect-ratio, SP current sheets. 
These were found to be violently unstable to the for- 
mation of plasmoid chains, the fastest growing wave 
number scaling as A: max Lcs ~ S 3 ^ 8 , with correspond- 
ing growth rate TmaxTA ~ g 1 / 4 , where Lcs is the length 
of the current layer, ta = Lqs/Va is the Alfven time 
(Va is the Alfven speed) and S is the Lundquist num- 
ber, S = LqsVa/v? where rj is the magnetic diffusivity. 
Since S ^> 1 in most applications of interest, this theory 
predicts the formation of multidinous plasmoids grow- 
ing super- Alfvcnically; the immediate implication is that 
stable reconnecting current sheets at large values of the 
Lundquist number cannot exist. These results have since 
been confirmed in direct numerical simulations [43l . l77j , 
and extended to account for the effect of a finite compo- 
nent of the magnetic field perpendicular to the reconnec- 
tion plane [78| and into the two- fluid regime (79j . 

The analysis of Paper I considered a very simpli- 
fied model background equilibrium, intended to retain 
only what we viewed as the most important features 
of a SP current sheet: a reconnecting magnetic field, 
B eq = (0,B y (x)) (x is the inflow direction, y the out- 
flow direction) and an incompressible flow defined by the 
stream function cf> eq — Toxy, where Tq = Va/Lqs is the 
flow shearing rate. The analytical derivation in Paper 
I did not, therefore, take into account potentially im- 
portant effects, such as the variation of the reconnecting 
magnetic field and of the outflow speed along the layer 
(i.e., along the y-direction in our chosen geometry), or 
the reconnected magnetic field. 

In this work, we generalize the results of Paper I to 
a more realistic, two-dimensional model of the current 
sheet. Using an approach in the spirit of WKB the- 
ory (justified by the expectation that the most unstable 
wave-number will be very large, fc max Lcs ^ 1)> we de- 
rive the dispersion relation for the plasmoid instability 
as a slow function of the position along the sheet, yo- We 
find that the scalings of the maximum growth rate (7max) 
and wave-number (fc max ) derived in Paper I hold true in 
a central, finite-sized patch of the current sheet; however, 
the growth rate and wave number are now parametrized 
nontrivially by yo- Surprisingly, we also discover that 
for a generic background equilibrium configuration, the 
maximum growth rate and wave-number of the instabil- 
ity increase with yo (i.e., outwards). As we show in this 
paper, a special point exists, j/o,criti beyond which the as- 
sumptions invoked in our calculation break down. This 
is the Alfven Mach point of the system, where the magni- 
tude of the outflow velocity (an increasing function of yo) 
becomes equal to the value of the Alfven speed based on 



the upstream magnetic field (a decreasing function of yo). 
Beyond that point, the current sheet becomes unstable 
to a different mode: the Kelvin-Hclmholtz (KH) insta- 
bility, whose growth rate and wave-number dependence 
we also derive analytically. 

The other main result of this paper is the study of the 
effect of a large viscosity v (parametrized by the magnetic 
Prandtl number Pm = v/r}^> 1) on the plasmoid insta- 
bility. The largc-Prandtl-number regime is pertinent to 
various astrophysical applications, e.g. , the interstellar 
medium (80j . and to fusion plasmas |3q |. and so it is im- 
portant to understand how large Pm affects plasmoid 
formation and dynamics. 

Our analytical results are complemented with a direct 
numerical solution of the full set of linearized equations. 

This paper is organized as follows. In section [Til wc 
present a heuristic derivation of our main results. A 
more rigorous approach to the problem begins in sec- 
tion IIII1 where the equations to be solved are laid out 
and the expected properties of the background equilib- 
rium are discussed (a more quantitative discussion of the 
constraints that the background equilibrium should sat- 
isfy can be found in Appendix [X[ where an analytical 2D 
SP-like current-sheet equilibrium is obtained). The core 
of the analytical calculation is presented in section IIVI 
The KH instability of the current sheet is derived in sec- 
tion |Vj Results of the direct numerical solution of the 
linear equations are presented in section IVII The effect 
of viscosity on the instability is addressed in section [VI Bl 
Finally, a discussion of the results and conclusions can be 
found in section [VTll 



II. HEURISTIC DERIVATION 

In this section, we show how the main results of this 
paper can be derived in a simple (albeit non-rigorous) 
way. A reader uninterested in the formal mathematical 
details can skip to section [Vl] after this section. 



A. Plasmoid instability 

The fastest growth rate of the plasmoid instability can 
be obtained from the usual tearing mode formulae as 

follows m. 

In the small A' limit, where A' is the usual tearing 
mode instability parameter, the standard (FKR) tearing 
mode dispersion relation is @ 
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(1) 



where a is the characteristic equilibrium magnetic field 
length scale, t v = a 2 jr\ is the resistive diffusion time and 
th = l/kE>o is the hydrodynamic time. In the opposite 
limit of large A' [8l|, 
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In both cases, the width of the resistive (or inner) bound- 
ary layer is 



^inncr/a ~ (l T H T r, *) 



(3) 



To find the fastest growing mode, let us assume a sim- 
ple Harris-sheet equilibrium, B y = Bq tanh(x/a); then, 
for fca <C 1, as we expect to be the case for k = fcmax, we 
have A'a ~ 1/fca. Substituting this expression in Eq. (TTJ), 
we find that it yields 7 oc fc~ 2 / 5 , whereas from Eq. @ we 
have 7 ~ A: 2 / 3 . Approximate expressions for the largest 
growth rate and corresponding wave number can there- 
fore be found by balancing Eqs. |T]) and (T5|). This gives 

k max a ~ (o/B ) 1/4 r- 1/4 , (4) 
Tmax ~ (a/B )- 1/2 r- 1/2 . (5) 
The corresponding inner-layer width is 

<Wer/a ~ (aBo/ry)- 1 / 4 . (6) 

In order to apply these scalings to a SP current sheet, 
we rescale the equilibrium length scale a to the sheet 
thickness: 



LcsS-V*. 



(7) 



Noting that the plasma outflow speed in a SP current 
sheet is Va = -Bo, we obtain 



Omner/oCS ~ S" 1 ^ 



(8) 
(9) 
(10) 



where ta = These predictions are in agreement 

with the results of Paper I [76[ • 

We can now use these results to estimate the critical 
value of the Lundquist number, S CI it, below which we 
expect SP current sheets to be stable. The underlying 
reasoning is that for the plasmoid instability to be trig- 
gered we must have 7 max T J 4 > 1, fc max L C s > 1 and 
<5inner/ocs <C 1, i.e., the instability has to grow faster 
than the characteristic outflow time and it has to fit in- 
side the current sheet, both along (thus the restriction on 
fcmax) and across (thus the restriction on <5i nnor ). Of all 
these conditions, the most stringent is that on the width 
of the inner layer, since it bears the weakest S depen- 
dence. Therefore, if we require (non-rigorously!) that 
<^inncr/^cs ~ 1/3 at the very most, then Eq. (fit)]) would 
yield SWit ~ 10 4 . This is consistent with numerical sim- 
ulations SIM El El. 



B. Plasmoid instability at large Pm 

One limitation of Paper I was that plasma viscosity 
was neglected. At low values of the magnetic Prandtl 
number, Pm = v /rj [relevant to the interiors of stars and 
planets, or liquid metal laboratory dynamos, for example 



[see [82| and references therein)] , the presence of viscosity 
should not change our results substantially. In contrast, 
for Pm S> 1, as is often found in fusion plasmas j3r| . 
warm interstellar and intracluster media |80l . [83j . etc., 
both the SP scalings and the tearing and kink modes 
scalings change [HJ|, and so, therefore, will the plasmoid 
instability. 

Let us work out the plasmoid scalings in the large- 
Pm limit in a similar way to that just presented for 
the invisicid case. Instead of the FKR @ and Coppi et 
al. [8l| results, we now use the corresponding formulae 
valid for Pm ^> 1, i.e., the the so-called "visco-tearing" 
(low A') and "visco-resistive kink" (large A') derived by 
Porcelli jH|. 

At low A' (the visco-tearing mode), we have 



7 ~ V 7 V'A'a, 

where r„ = a 2 /v is the viscous diffusion time. At A' 
oo (the visco-resistive kink), the growth rate is 
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_ -2/3 -2/3 1/3 
I H T ri T v ■ 

The corresponding inner-layer width is 



Omncr/a ~ [iff / (TjjT,,)] 



1/6 



(12) 



(13) 



As before, let us assume that A'a ~ 1/fca for k ~ 
fcmax- Then we find from Eq. (pTT|) that 7 oc fc~ 2 / 3 , 
whereas Eq. (fl"2")) yields 7 oc fc 2 / 3 . Scalings for the fastest 
growing mode can thus again be found by balancing Eqs. 
(fnj) and (fTS|>. The result is 

fc m a X a~(a/So) 1/4 T- 1 / 8 r- 1 / 8 , (14) 

7 max -(a/So)- 1/2 r- 3 / 4 ry 4 , (15) 

<W/a ~ (a/B ) 1/4 r- 1/8 r- 1/8 . (16) 

We now repeat the previous procedure of rescaling the 
equilibrium length scale a, this time using the results 
obtained by Park et al. [36| for the SP model in the limit 
Pm > 1: 



This gives 



a = ocs ~ LcnS-^Pm 1 / 4 . 



fcmaxicS ~ S 3/8 ^" 3/16 , 
7ma X ^^^/ 4 Pm- 5 / 8 , 

<5in„cr/<5cs ~ S-^Pm 1 ' 16 . 



(17) 

(18) 
(19) 
(20) 



We shall find in section I VII that these scalings indeed 
agree very well with the results of a direct numerical inte- 
gration of the linearized equations. We thus find that the 
dependence of 7 ma x and fc max on S remains unchanged 
at large Pm. However, viscosity damps the instability, 
and decreases the wave-number and the growth rate of 
the fastest growing mode, while slightly thickening the 
inner layer. 
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An important question is how the critical Lundquist 
number for the onset of the current sheet instability, 5 cr ; t , 
scales with the magnetic Prandtl number. Although the 
expressions above are formally only valid in the limit 
S ^> 5 cr itj P m 3> 1, we can use them to obtain a 
rough estimate of this dependence. Since the instabil- 
ity requires <5inner/^cs *C 1, we may again demand that 
this be at most 1/3 and use Eq. ([2T>1) to obtain 



10 4 Pm 1/2 . 



(21) 



To see the consistency of this result, note that the 
same dependence of S ct it on Pm can be obtained by 
either looking for the minimum wavenumber that will 
fit inside the current sheet, fcmax-^cs ~ 1, or by re- 
quiring that the growth rate is comparable to the flow 
shear rate, 7 max £csAwt ~ 1 (note that for Pm 1, 
Wout ~ V^Pm -1 / 2 |36j|). In both cases, Eq. ((TSj) and 
Eq. (|T^|) yield S C rit ~ Pm 1 / 2 . This result is a specific 
prediction, which can in principle be checked via direct 
numerical simulations of the current-sheet instability in 
the large-Prandtl-number regime. 



C. Kelvin- Helmholtz instability 

The velocity outflow profile of a Sweet-Parker recon- 
nection configuration is such that it is maximum at the 
midplanc (x = 0) of the current sheet, and decays to zero 
away from it. Thus, there are two parallel shear layers, 
with two corresponding inflection points of the outflow, 
at x ~ ±<5cs- The shear (in the ^-direction) of this flow 
profile can be estimated as: 



du,. 



v A y 



dx 5cs Lcs 



(22) 



Each of these layers would be Kelvin-Helmholtz unsta- 
ble were it not for the stabilizing effect of the upstream 
magnetic field, B y (85|: as is well known, a magnetic 
field that is coplanar with the flow profile will stabilize 
the KH instability as long as \B y \ > \u y \. In the case 
of a SP current sheet, the upstream magnetic field B y is 
not constant along the sheet; in particular, its magnitude 
decreases in the y-direction; a simple model for it is [86| 



B„ — B, 



2 

cs- 



(23) 



It is thus possible that there exists a location j/crit/^cs ~ 
1 along the sheet where the (decreasing) strength of 
By matches the (increasing) magnitude of the ouflow, 
Uy ~ Va y/Lcs- This is the Alfven Mach point of the 
system; for y > y C mt, the magnetic field is no longer able 
to stabilize the KH mode. Fig. [1] provides a schematic 
illustration of both types (KH stable and unstable) of 
configuration. 

A rigorous derivation of this instability is presented 
in section [V] Here we show how the basic scalings for the 
fastest growing mode and corresponding wave number 




FIG. 1. Cartoon illustrating KH-stable and KH-unstable 
parts of an idealized SP current sheet. The outflow profile 
u y is depicted in blue; the upstream magnetic field B y in red. 
The vertical dashed lines mark the position of KH-unstable 
layers (i.e., the inflection points of the outflow profile, where 
d 2 u y /dx 2 = 0). The full lines show the profiles below the 
Alfven Mach point (y < y cr it); in this case the magnitude 
of By at x ~ ±5cs exceeds that of u y at x — 0, and thus 
the magnetic field is sufficiently strong to stabilise the cur- 
rent sheet against the KH instability. The opposite case of 
profiles above the Alfven Mach point ( y > y C rit) is shown by 
the dashed lines; B y at the inflow is now weaker than u y at 
the center of the sheet, which is, therefore, KH-unstable. 



can be obtained heuristically. From the standard theory 
of the KH instability [8f| , it is known that 



kh du v 
7raax ~ dx 



Y± 1.KH 



1. 



(24) 



where a is the characteristic scale length of the sheared 
flow profile u y (x). As above, these estimates can be 
applied to the SP-sheet by simply setting a — > 6qs ~ 
LcsS" 1 / 2 , implying that fc™ - l/tfcs- Thus, 



7 KH r, 

/ max 



L 



cs 



(25) 



Note that this growth rate is even larger (i.e., has a 
steeper positive scaling with S) than that of the plas- 
moid instability, Eq. ©. 

The KH instability scalings in the large Pm limit 
are easily obtained in a similar way, using the follow- 
ing modifications of the SP relations derived in [3t| : 
Mout - VaPtu- 1 ' 2 , S cs ~ LcsS^/iPm 1 / 4 . Then 
we have 



I max & 



5Va Pm -s/4 



k KH Lr« 



(26) 



III. PROBLEM SETUP 



In this section we proceed to make some of the above 
discussion more rigorous and quantitative. 
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We solve the 2D rcduccd-MHD equations [87 



{V,ViV}+^Vl0, (27) 



(28) 



Here <j> and V are the stream and flux functions of 
the in-plane velocity and magnetic field, respectively, so 
u = (—d y (j), d x 4>)i B = (—dytf), d x ip); the magnetic field is 
measured in velocity units; Poisson brackets are denoted 
by {</>, tp} = d x 4>dyip—dy(f)dxtp; T] and v denote the plasma 
resistivity and viscosity, respectively; Eq represents an 
externally applied electric field, required to sustain an 
equilibrium in the presence of finite resistivity. In the 
analytical calculation that follows, we will assume that 
the magnetic Prandtl number is small, Pm = vjr\ <C 1 
and therefore neglect the effect of viscosity on the linear 
instability. The case of large Pm will be studied numer- 
ically in section fVI B[ 

We are interested in analysing the linear stability of 
an SP-like current sheet, whose inverse aspect ratio, e, is 
predicted by the SP model to scale as 



e = 8 C s/L 



cs 



-1/2 



< 1. 



(29) 



Therefore, in the vicinity of a general point along the cur- 
rent sheet, y — yo, and provided that yo is not too close 
to either of the ends of the sheet, it is reasonable to ex- 
pand the equilibrium magnetic flux and stream functions 
in a power series (88j : 



4>eq (X,y) \y=y = ^ — } 



77=0 



"I 



(X,y)\y=y = 



77=0 




ipn(x,yo), (30) 



<f>n(x,Vo). (31) 



The functions tpn{x,yo)i (pn( x ,yo) can in principle be 
found either by substituting these expansions back into 
Eqs. (f2"7H2"5|) and solving the equilibrium problem order 
by order in (y— yoV-^cSi or by Taylor expanding a known 
equilibrium solution, ip eg , <f) eq around the point y = yo. 
Neither of these procedures is straighforward: the first 
option implies either truncating the expansion at some 
arbitrary order or guessing one of the equilibrium func- 
tions [88| to solve the closure problem; the second re- 
quires an exact analytical solution. To the best of our 
knowledge, a 2D analytical solution to this problem that 
would capture all the essential features of a resistive SP- 
like current sheet equilibrium has never been derived [8!| . 
Fortunately, however, we shall find in the following sec- 
tions that, actually deriving the linear instability requires 
very little information about the equilibrium profiles, and 
the problem can be solved for general functions ip n , <f) n 
provided that the following key assumptions hold: 

• the background equilibrium observes the expected 
symmetries, i.e., at the center of the sheet 

V>2»+i(2/o = 0) = 0, </> 2n (y = 0) = 0; (32) 



• for |x| 3> 1, the incoming flow u x and the recon- 
necting magnetic field B y approach constant (in x) 
values (which can, however, be functions of yo)- 

Further general properties of the equilibrium that will be 
needed in our calculation are derived in Appendix |A"1 

Note that, in Paper I, we considered only the case 
yo = 0. Furthermore, we adopted a very simplified 
description of an equilibrium current sheet in which only 
<j>i and tpo were non-zero, with <f>i = Toxy, tpo = 
i.e., the flow had no vorticity ui z = V^^i = 0, and 
the reconnected field was ignored, B x = —dijj/dy = 0. 
In this paper, we drop those model assumptions and 
consider an arbitrary, two-dimensional current-sheet 
equilibrium. 



A. Normalizations 

We introduce the following normalizations, motivated 
by the SP scalings: 



<t>'o = r oyow(x); 

4>i = -T Q LcsScsu(x); 

<f>2 = ~r ?/ <5csM>(x); 



= B f(x); 
i'i = -T y S C sg(x); (33) 
i'2 = -T L C sScsh(x); 



where 



£cs = {ri/To) 1/2 



(34) 



and To = Bq/Lqs (note that in our units Bo = Va)- Wc 
also normalise time and lengths as follows: 



tr = t; x/<5cs=£; y/L C s = y- 



(35) 



Under these normalizations, the magnetic and velocity 
fields obtained from the power-series expansions (j30H3Tj) 
keeping only up to first-order corrections in y — yo are: 

B y I Bo = /(C) ~(y- yo)y~og'(0 (reconnecting), (36) 
B x /B = e [y g(0 + {y~ yo) HO] (reconnected) (37) 



and 



u y /V A = y v(0 -(y- yo) "'(£) (outflow), (38) 
u x /V A = e [u(0 - (y - y ) y Q w(0] (inflow). (39) 

It is clear from these expressions what the physical sig- 
nificance of the functions /, g, h, u, it, w is. The phys- 
ical units and presumed magnitudes of these fields have 
been absorbed into the normalizations, so these functions 
are all order- unity dimensionless quantities. They can 
have parametric dependence on yo, but note that the pre- 
sumed lowest-order linear dependence of the reconnected 
magnetic field and of the outflow on yo are explicitly in- 
cluded in the normalizations Eq. (|33[) . 
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B. Linearized Equations 



Let us consider small perturbations to a generic equi- 
librium, ip = ip eq + 6tp(x,y,t), 4> = (f> eq + 6<j)(x,y,t), 
and linearize the RMHD Eqs. (|27ILT51) using the expan- 
sions (|30H3ip for the equilibrium profiles, keeping terms 
up to first order in (y — y ). We obtain 



^ + [yov(0 -($- y ) «'«)] ^ + HO + {v- wo) vMO\ 
- [yog(0 + (y- m) HO] ^ - [/(f) - (y - yo) yog'®} ^ = 



dy 2 



# + o((y-y ) 2 ), (40) 



Y T + [5b«(0 - (y - w>) ^ + [«(f) + - yo) y<Mf)] ^} + ^ 

- [«"(0 + (y - yo) f/oV'(0] ^ - [vov"(0 - (y yo) </"(£)] |? 
[f(0 - (y - m) yog'm + [yog(0 + (y - yo) H0\ ^} + fi2 ^) W 



&7bs"(0 + (y - So) h"(0] ^ - [/"(£) + (y - y ) yog"] ^ + O ((* - yo) 2 ) 



9f 



<9y 



(41) 



where we have used the normalizations defined in 
Eqs. d33H35]). 

In the case of yo = 0, the above equations include the 
following effects which were absent in Paper I: 

• in Eq. ([4TJ]). the term proportional to h(£), repre- 
sents the effect of the reconnected magnetic field; 

• in Eq. (|4"Tj) , the terms proportional to u"(f) and 
u"'(£) on the left-hand side, which represent the 
vorticity of the equilibrium flow; 

• the term proportional to h"(£) on the right-hand 
side of Eq. (|4"Tj) , which is the contribution to the 
equilibrium-current gradient from the reconnected 
magnetic field. 

Further progress at this point is hindered by the fact 
that these equations contain explicit dependences on the 
y variable and cannot, therefore, be Fourier transformed 
in this direction. To address this difficulty, let us compare 
the magnitudes of the first and third terms on the left- 
hand sides of these equations: 



d/dr 
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(y-yo) d/dy (y-yo)^' 



(42) 



where 7 is the growth rate of the anticipated instability 
at y = yo normalized to the Alfvenic shearing rate To 
and 

K = kL cs (43) 
is the normalized wavenumber of the perturbation at that 
location. Thus, the third term can be ignored if the anal- 
ysis is restricted to patches of the current sheet whose 
extent in the y-direction is such that 



(y - yo) < i/k. 

This approach is valid provided that 

K (y - yo) > l- 



(44) 



(45) 



In other words, the domain in the y-direction is divided 
into smaller patches, and the linear analysis performed 
locally in of each of these. A WKB approach remains 
valid provided that asymptotically many wave lengths fit 
in each of these patches. Eqs. (|44ft45|) imply that we seek 
solutions such that 



d_ 



» 1, 



d_ 

dy 



> 1. 



(46) 



These are a priori assumptions, which will be later jus- 
tified by our ability to obtain such solutions. 

Under these approximations, Eqs. (|4"0H4"Tj) become 



Note that the functions h(£) and u>(£) have dropped out. 

One can now look for linear modes of the form 
exp(— iujr). In the midplane of the current sheet (i.e., 
£ = 0), the plasma is flowing outwards at some fraction 
of the Alfven speed, yovo, where vq = u(£ = 0, yo). Let 
us take this into consideration explicitly and set 



ui = ny v + i-f. 



(49) 



We therefore look for solutions to Eqs. (|4"TH4"5|) in the 
form: 



Sip = *(£)e 
5<j> = 



[■y-iKy v ]T+iKy 



,\y—iKy vo\T+iK,y 



(50) 
(51) 



We are ignoring the time dependence of k due to the 
background flows, because this variation will occur on a 
much longer timescale than that of the expected growth 
rate of the instability |h| • Introducing the parameter 



A = 7/k, 

we obtain 

{a-^ [« -«(£)]}* + ^ 



(52) 



(53) 



{A -iy [«b -«(€)]}(*" 



K 2 e 2 $) 



K k 
— /(O (*" - * 2 e 2 *) + i^-g{Z){y"' - k 2 eH') 



i^ g "m'+f"m- 



(54) 



This is the set of equations that will be solved in sec- 
tions EH and El 



IV. PLASMOID INSTABILITY 



We will now proceed to solve Eqs. (|5jH54)) in three 
different regions: the "external" (global) region where 



<5 C s < % < L (i.e., 1 < £ < S^ 2 ), the "outer" region 
(the SP current sheet), where x ~ <5cs (£ ~ 1) and fi- 
nally the "inner" region (the inner layer inside the current 
sheet where plasmoids form), where x <C <5cs (£ *C 1). 
Outside the current sheet (i.e., in the external region), 
the plasma is ideal; inside the current sheet (i.e., in the 
outer and inner regions), resistive effects cannot be ne- 
glected. In addition to the assumptions of Eq. (|4*5)l . we 
will require here that A = 7/ ft <C 1. This ordering is 
indeed satisfied by the fastest growing mode, as we will 
confirm a posteriori. 



A. External Region: |£| » 1 

This is the upstream region outside the current layer, 
i.e., x ^> 5qs- Here, we expect the equilibrium profiles to 
behave as: 



«(0 -> 0, 

/(0 « ±/o 



«(0 

5(0 



(55) 
(56) 



where 

^00, /cxd , 5oo are functions of yQ only, taken to be 
of order unity (see Appcndix|X]); of these, we will discover 
that only matters for the calculation of the instability. 
In the above expressions, the upper sign applies to £ > 0, 
and the lower sign to £ < (so as to observe the expected 
parities of the equilibrium, namely that it(£) and /(£) are 
odd in £, and <?(£) is even). The linear dependence of g(£) 
on £ for large £ might not be obvious at first glance and 
is derived in Appendix [AJ 

We will make the a priori assumption that the terms 
proportional to u(£) or and to <?(£) or <?"(£) in 

Eqs. ([53TI54"]) are negligible in this region, and then show 
that this is indeed the case. In the absence of these terms, 
we obtain: 

(A - H/0«o) * = ±/o«,*, (57) 
(A - iy v ) - ft 2 e 2 $) = T/oc(*" - « 2 e 2 *), (58) 

which can be easily combined to yield: 

[fl - (A - iy v Q ) 2 ] - ftV*) = 0. (59) 
The general solution to this equation is simply 

* ± = e*"* (60) 
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where C3 are integration constants and ± refers to £ ^ 
[in unsealed units: = exp(=Ffea;)]. 

We can now check the assumption about the smallncss 
of the terms proportional to <?(£) and <?"(£)■ 

From Eq. (|60[) and Eq. ([57]) we see that, for arbitrary yo, 
we have \?'/\& ~ $'/$ ~ Ke. Let us then compare the 
magnitudes of the third and first terms on the left-hand 
side of Eq. ([55)1 [the same reasoning applies to Eq. : 



A* 



A" 



(61) 



which is small provided that A 3> e, a condition we will 
later see is satisfied by the fastest growing mode. With 
respect to terms involving <?(£), the ratio of the magni- 
tudes of the fourth and fifth terms on the left-hand side 
of Eq. (f53"f [and similarly for the second and first terms 
on the right-hand side of Eq. ([53)1] is 



B. Outer Region: |£| ~ 1 

This region represents the SP current sheet itself, i.e., 
x ~ Sqs- Here, the functions /(£) and g(£) 

are simply assumed to be ~ 0(1). We find that terms 
proportional to u(£) and <?(£) or to their derivatives are 
again negligible. For example, consideration of the same 
terms as in the previous section leads to: 



- 7 1 < 1 



A* 

yog® I* 



y K 1 < 1, 



(63) 
(64) 



and similarly for the others. Therefore, to lowest order 
in e, Eqs. (J5H] become 



Wo 



(62) 



This is again small provided that £ <C e _1 ~ we 
shall find that the fastest growing wave number K max ~ 
e~ 3 / 4 and thus the eigenfunction, Eq. decays before 
the condition £e -c 1 breaks down. Note also that the 
expression for <j(£) given in Eq. (|56]l is not expected to 
hold for £ > e _1 (i.e., x > Lcs)- the reconnected mag- 
netic field does not grow unbounded as £ — > 00). 



{\-iyo [«b -««)]} * = /(£)*> 



{A - ^y [w - «(£)]} (*" - K 2 e 2 $) - iy v"(0® 

+ [no + ^ 2 /(o] 



(65) 



(66) 



Combining these equations results in the following eigen- 
value problem: 



no 

/(0 



2 2 

« e 



A - iyo bo - v(Q] 

/(o 



de 



2 2 

k e 



A - iyo [v - v(£)} iyov"{0 



/(0 



/(0 



(67) 



r 



subject to boundary conditions given by the external so- 
lution, Eq. (|60p . and the requirement (for the plasmoid 
instability) that \1/ be an even function. Writing the so- 
lution in the form = /(0x(0) Eq. (|6"?| becomes: 

^{V(0x'm=t 2 K 2 V(0x(0, (68) 
where the "potential" V(£) is 

nO = / 2 (0 + [A - i yo(«o - «(£))] 2 ■ (69) 

Note that in the case yo -C 1 all terms on the right- 
hand side of Eq. (|67| are small compared to the first 
two terms on the left-hand side. In that case, we recover 
to the problem solved in Paper I [except here /(£) re- 
mains unspecified]. For the general case yo ~ 1, an exact 
solution can be obtained provided that the terms propor- 
tional to K 2 e 2 can be neglected (we will later check the 
validity of this assumption; for now, let us call attention 
to what it means: fte = kdes, so the assumption ne -C 1 
implies that the wavelength of the expected instability 
is much longer than the current sheet thickness). Thus, 



we neglect the right-hand side of Eq. (f6"8")l and find the 
solution 



no' 



(70) 



where ± refers to £ ^ 0, C^, are constants in in- 
tegration and £0 is an arbitrary number of order unity 
(different choices of £0 will produce subdominant correc- 
tions to Cf'). 

For the plasmoid instability, we expect A < 1, so this 
solution simplifies to 



(71) 



We now match this expression to the external solution, 
Eq. (|60)) . in the region 1 <C £ <C (ne)^ 1 , or, equivalently, 
in dimensional form, 5c s *^ x <C 1/k. In this region, 
«(±0 « v and f'(±0 < 1, implying that |/(±£)| a 
and the integral is dominated by the upper limit. We 
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obtain 



r ± _ _ 



C*2 /oo 



-2 2 ' 



(72) 
(73) 



For £ « 1, /(C) « /5C and = v - <£ 2 /2. Thus, 
to lowest order in £, the integrand in Eq. (|71l) becomes 
1/ (/q £ 2 ) . The integral is therefore again dominated by 
the upper limit and we obtain: 



*±(0) = -£_. 

Jo 

Demanding that be an even function of £ yields 

ci = c 2 - = -/^>(0). 

Thus, the outer region solution is: 

1 1 



(74) 



(75) 



= /o*(0)/(0 ± 



" /2 

J OO 



y 2 u 2 Ke 



(76) 



Note that the second term of the above expression en- 
sures that the eigenfunction remains finite at £ = 0. 

As usual in tearing-mode-type calculations, let us now 
introduce the standard instability parameter 



A' 



*'(+0) - *'(-0) 



*(0) 

Using Eq. ([76]), we obtain (§3 



A' 



Jo 



-2 2 ' 



(77) 



(78) 



We stress that the functional dependence of A' on yo 
is both explicit and implicit, as /q, /oo and vq are all 
functions of yo. 



C. Inner Region: |£| < 1 

This region is the internal layer inside the SP current 
sheet, i.e., x <C Sqs- Wc begin by noting that, again, 
independently of the specific functional form of the SP 
equilibrium, the symmetries of the problem are such that, 
for |£| <C 1, the equilibrium profiles can be approximated 
as 

f(o = fot+o(e), g(o=.9o+o(a, (79) 

u(0=^ + O(C 3 ), v(Z)=v + O(e), (80) 

where /q, u' , go, vq are constants with respect to £ but 
depend on yo- In this region, the relative magnitudes of 



the different terms in Eqs. 
of the following cases: 

u*' 1 

K A* 

yog®' i 

K /$ 
K A* 

u"$' 1 



can be reduced to one 



Whinner 'I' K 



7* 



1 

■ 

7 



1 

7^ 

5: 



1X2 

inner 



k8?. 



inner 

2 

inner 



K A$" 

A$" 

yog"*' i 
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1 

K 

Except for Eq. (|83| . all these ratios can be shown 
teriori to be small. Thus, to lowest order, Eqs 
become 

a* - yfc* = 



(81) 

(82) 

(83) 

(84) 

(85) 

(86) 
a pos- 



A$' 



-/of*"- 



(87) 
(88) 



These equations are mathematically the same as the 
equations for the tearing mode in the inner region, ex- 
cept that here the role of the small parameter is played by 
1/k rather than resistivity. Since A'<5i nnor is not expected 
to be small, the constant-* approximation Q cannot be 
used. Instead, this eigenvalue problem is mathematically 
equivalent to the one solved by Coppi et al. for the re- 
sistive internal kink mode [8lj . The resulting dispersion 
relation is 



4(*/o) 1/3 A 5/4 



r [(a 3 / 2 - 1) /4] 

r[(A 3 / 2 + 5) /4] 



= A' 



Jo 



(89) 



-2 2 ' 

Vo v o 



where T is the gamma function and A = 7/g 2 ^ 3 n 2 / 3 . 
The width of the inner region is 



= ^1/^-1/2^-1/2 



7 ' K 



(90) 



which can be used to confirm the smallness of the terms 
neglected in deriving Eqs. (|87H55]) . 

D. Solution of the Dispersion Relation, Eq. (|89p 



The dispersion relation, Eq. (|89f , has two relevant lim- 
its. For A < 1, 



0.95- 



16 T(5/4) 

Tr(-i/4) 
fi - 

Jo 



4/5 



Jo 



2„,2\ 4 /5 



2„,2\4/5 



K - 2 / 5 e- 4 / 5 . 



(91) 
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This is valid provided that k 3> e~ 3 / 4 (but also ne <C 1, 
as required by our earlier assumptions). 

On the other hand, taking A — > 1— (from below), we 
obtain 

7 = (/^) 2/3 -^ / ~"fr 2 « 2 e , 02) 
J Jo 

valid for k < e~ 3 / 4 . The scaling of the fastest growing 
wave number can be determined by balancing the two 
terms on the right-hand side of the above expression. 
The exact result can be obtained by solving the equation 
dy/dn = in the limit A — > 1— . We obtain 

( 1 /o 5/3 \ ^ -3 4 

Kmax = [v^Tl^Wo) e_3/4 ' (93) 

7maX= 3^V/--^o 2 ^ 1/2 ' (94) 

The scalings with e are the same as those derived in 
Paper I, although here they have been obtained for a 
general SP equilibrium. In other words, this linear the- 
ory predicts that current sheets are unstable to a super- 
Alfvenic instability whose growth rate increases with the 
Lundquist number, 7 max ~ S* 1 / 4 [recall that e = S* -1 / 2 ; 
see Eq. (f2T)|) ]. A plasmoid chain forms inside a region of 
width dinner ~ S~ 1 / 8 6qs, with the number of plasmoids 
scaling as K max ~ 5 3 / 8 . These scalings justify the order- 
ing assumptions employed in deriving these results, i.e., 

7max > 1, Kmax > 1 and A 

max ^ 7max / ^max ^ 6 ^ , SO, 

e < A max < 1. 

Let us now analyze the dependence of 7 max and K max 
on the position yo along the sheet. Note first that the in- 
stability vanishes at the locations where the equilibrium 
current Jq = — the end points of the SP current sheet. 
Note also that since we have dropped corrections of or- 
der e 1 / 2 in our derivation [for example, the terms propor- 
tional to K 2 e 2 in Eq. the terms proportional to y$ 
in the above expressions are only to be kept if yo 3> e 1 / 4 . 
For values of yo ^ e 1 ^ 4 , all yo corrections are negligible 
to lowest order and Eqs. (|9"3TIM)) simplify to yield 

W = vr- 3 / 8 £ 5 /V 3 / 4 , (95) 
7max = ^ £ 3/ W 2 , (96) 

where Eq = LcsEq/(B 2 Scs) is the normalized back- 
ground electric field, and we have used the relation- 
ship /q = Eq, which follows from Eq. (|A3|) in the limit 
Vo < e 1/4 [note that f^yo = 0) = 1]. 

In the opposite case of yo ^> e 1 / 4 , the dependence of 
K max and 7 max on yo is a nontrivial function of the spe- 
cific values of the equilibrium coefficients, which are all 
functions of yo (this does not affect the scaling of « ma x 
and 7 max with e). For y <C 1, exact values of the coefi- 
cients of the Taylor expansion of the equilibrium around 
£ = were derived semi-analytically by Uzdcnsky and 



Kulsrud [92j assuming a Syrovatskii-like upstream mag- 
netic field [Eq. (|23]) ] see Appendix |A"1 Using those co- 
efficients and the relationship /q = Eq — ypWoffo derived 
in Appendix [3] a Taylor expansion of Eqs. (j9"3TIM|) in yo 
yields: 

>W(N « 1) « (0.56 + 0.10y 2 ) e~ 3/ \ (97) 
7max(|yo| « 1) « (0.71 + 0.73y 2 ) e^ 2 . (98) 

These results reveal a perhaps unexpected feature of 
the plasmoid instability: both K max and 7 ma x increase 
with distance from the center of the sheet. The same 
conclusion is easily deduced for arbitrary yo ~ 1: even 
though in that case one is forced to retain the full ex- 
pressions for K max and 7 max , Eqs. (|9"3"HM)) . it is also 
true that, under very general conditions, one expects /oo 
to be a decreasing function of yo (for example, a stan- 
dard Syrovatskii ideal current sheet solution [86[ yields 
foo = \/l - Vq; see Appendix EJ). 

A problem thus arises: it is possible that a loca- 
tion yo = yo,crit ~ C(l) exists inside the current sheet 
where /^ — y~Q crit i>o = 0, and our solution breaks down 
(A', 7, k — > +oo). It is clear that while approaching 
that point, K max will get to be so large that terms of 
order K max e, or Kniax e2 > can 110 longer be neglected and 
thus our ordering assumptions become invalid. For val- 
ues of | j/o | > yo.crit; our solution is again physical, but 
the second term on the right-hand side of Eq. (|92|) equa- 
tion changes sign, implying that 7 grows with k and the 
value of the fastest-growing wave number cannot be de- 
duced from this equation. This means that terms propor- 
tional to K 2 e 2 are still necessary to determine the fastest 
growing mode. Physically, the fact that the condition 
K max e ~ 1 is met somewhere in the sheet means that at 
those locations the wavelength of the perturbation be- 
comes comparable to the current sheet thickess, whereas 
for 1 2/0 1 "C yo,crit it was much longer. 

It is easy to understand why the location — 
2/ocrit^o = should be special: this is where the mid- 
plane outflow speed (i.e., u y measured at x = 0) matches 
the Alfven velocity associated with the upstream mag- 
netic field (measured at x ~ <5cs); i.e., it is the Alfven 
Mach point of the system. The background outflow ve- 
locity profile is strongly sheared (in the x-direction) in- 
side the current sheet. Such a profile would be unstable 
to the Kelvin-Helmholtz (KH) instability, were it not for 
the stabilizing effect provided by the (flow-aligned) back- 
ground magnetic field B y (85| . However, whereas the flow 
grows in magnitude with increasing y, B y does the oppo- 
site (see the discussion of section lll C|) . The Alfven Mach 
point is where the two amplitudes match. Beyond that 
point, the magnetic field is no longer sufficiently strong to 
provide stability, and we should thus expect the plasmoid 
instability to morph into the KH instability. The increase 
of 7max and K max along the sheet is thus a reflection of 
the fact that the current sheet is increasingly unstable to 
the KH mode. For yo > yo.crit the plasmoid instability 
is replaced by the KH instability as the most unstable 
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mode. The most unstable wave number also grows along 
the sheet because, as is well known [85[ , the growth rate 
of the KH instability peaks at ne ~ 1 (i.e., kScs ~ lj 
where 5cs is the scale of the cross-sheet velocity shear). 

To summarize, the analytical dispersion rela- 
tion, Eq. accurately describes the plasmoid 
instability of the current sheet for — j/o,crit < Vo < 2/o,crit- 
For values of y~Q outside this interval, it becomes nec- 
essary to keep terms proportional to n 2 e 2 in order to 
calculate the fastest growing mode, which is no longer 
the plasmoid instability, but the KH instability. Re- 
taining the n 2 e 2 terms analytically is difficult; however, 
we have been able to address the full problem by a 
direct numerical solution of the linearised equations. 
Results are shown in section fVTl Before discussing those, 
though, let us present an analytical derivation of the KH 
instability of a current sheet valid in the long- wavelength 
limit ke< 1. 



V. KELVIN-HELMHOLTZ INSTABILITY OF 
THE LAYER 

In this section, we derive an analytical dispersion re- 
lation of the Kelvin-Hclmholtz (KH) instability of the 
current sheet valid for perturbations whose wave number 
is k <^ e _1 . This ordering of K is not expected to capture 
the fastest growing mode, as suggested by the heuristic 
derivation of section Til CI (K max ~ e _1 ), but we are un- 
able to obtain an analytic solution valid for ne ~ 1 . The 
derivation presented here, however, does reveal a number 
of interesting features of the KH instability of the cur- 
rent sheet. The direct numerical solution presented in 
the next section does of course cover all values of k. 

We first remind the reader that in a SP reconnection 
configuration, the outflow velocity profile is maximum at 
the midplane (x = 0) of the sheet and decays to zero 
away from it. There are, therefore, two shear layers on 
each side of the sheet, at x ~ ±<5cs> where the KH insta- 
bility may develop — see Fig. [T] These two shear layers 
will push magnetic fields of opposite sign on each side 
of x = towards each other. This will create a current 
sheet at x = 0. Thus, resistivity can play an impor- 
tant role in this mode, not at the KH layers themselves, 
x ~ ±<5cS; but at x = 0, where the shear layers interact. 
This situation is then conceptually similar to the forced 
reconnection problem treated by Hahm & Kushud (93j 
(the Taylor problem), where perturbations at some dis- 
tant boundaries on each side of a rational surface induce 
reconnection at that surface. Here, the KH instability 
can be thought of as the equivalent of those perturba- 
tions at the boundaries, forcing reconnection at x = 0. 

In [93| (see also (94[) it was found that if the change 
in the boundaries occurs on a timescale much faster than 
the resistive one, there will be no reconnection in the 
early (linear) stage of evolution; the magnetic field will 
pile up until the current gets sufficiently large for the re- 
sistive term to become important. In a somewhat similar 



fashion, here, whether or not the KH instability induces 
reconnection depends on how large its growth rate is; for 
low values of ne, when its growth rate is lower, the pile 
up of the magnetic field is prevented by reconnection, 
which can proceed at a rate comparable to the growth 
rate of the KH instability. However, for larger values 
of Ke, the KH instability will be faster than the rate at 
which reconnection can occur and we expect to find an 
ideal (i.e., non-reconnecting) mode (though reconnection 
is expected to start in the nonlinear stage, which we do 
not address here). 

Our formal analysis of this problem again considers 
three asymptotic regions: the inner region (|£| <IC 1), the 
outer region (or flow shear layer, |£| ~ 1), and the exter- 
nal region (|£| ^> 1). Since all the same orderings that 
led to the simplification of Eqs. (|53HMf in the previous 
section are still expected to hold here, the equations to 
solve in each of these regions are the same. In partic- 
ular, the solution in the the external and outer regions 
remains unchanged and is given by Eqs. (|60p and (|7ip . 
respectively. 

In the inner region, wc must solve Eqs. (|57H55|) . There 
are two cases of interest: when the right-hand side of 
Eq. (|87[) is important, and when it is not. The former 
case occurs at low values of Ke, whereupon we simply re- 
cover the dispersion relation Eq. (|9"2"]l . This is similar to 
the resistive- kink-mode solution of Coppi et al. [8lj : as yo 
increases and eventually becomes such that yoVQ > f^, 
A' [Eq. ((75)) ] transitions from positive to negative via 
infinity; mathematically, this is equivalent to the well 
known transition from the very unstableQarge-A') tear- 



ing mode to the resistive kink mode |81| . except here 
reconnection is driven by the KH instability, rather than 
by the kink mode. 

For larger values of tie (though still requiring that 
K€ -C 1 ; we will later determine how large K has to be for 
the following to hold), we can ignore the right-hand side 
of Eq. dSJJ and obtain 



<lK 



So, Eq. 



becomes 



(99) 



(100) 



to be solved subject to the boundary condition $(0) 
(H. The solution is: 



$ ± (e) = ^arctan(4c), 



where $q = <1>'(0). The inner layer width is, therefore 

A 



(101) 



8h 



ft 



(102) 



Using Eq. (TToTj) to substitute for $ in Eq. ([99]). we 
obtain 



* ± (0 = ^arctan( ^ 



(103) 
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As expected this is a non- reconnecting mode: "5(0) = 
— indeed, this solution is mathematically equivalent 
to the ideal-kink-mode solution found by Rosenbluth et 
al. [§6| ; physically, the difference here is that the drive is 
the KH instability. 

The solution, Eq. (|103p . is now matched to the solu- 
tion in the outer region, Eq. ([7I]). For £ > 1, Eq. (fT05)) 
becomes: 



A 



Therefore, we have, 



Ct = A$' . 



(104) 

(105) 
(106) 



Substituting these expressions for Cf, C^ in Eq. ([75]). 
we obtain the final dispersion relation: 

x ^ yo v O ~ foe 2 / in7 \ 

7 = Ak=- k e. (107) 

Z JQ 

This expression shows that an unstable mode exists when 
\Vo v o\ > I/oo|j i-e., above the Alfvcn Mach point of the 
system, \y~o\ > Vo.mt- To determine its region of valid- 
ity, let us use Eq. (|107[) and Eq. (|102[) to compare the 
the right-hand side of Eq. ([57]) , which we neglect in this 
derivation, to the first term on the left-hand side of that 
equation. We then find that this solution is valid for 
6~ 3/i < k < e _1 . 

Finally, we estimate the value of n = n tr above which 
Eq. (|107j) yields faster growth than the resistive disper- 
sion relation, Eq. (p?2"]) . This corresponds to a transition 
from a reconnecting to a non-reconnecting mode. Com- 
paring the two expressions we find 



-3/4 



f 



5/4 



3/4 



-3/4 



(108) 



Wave numbers such that k > K tr are ideal, non- 
reconnecting modes (i.e., ^(0) = 0). This implies, in 
particular, that the fastest growing mode, K max ~ e , is 
an ideal mode, as will be confirmed by the full numerical 
solution presented in the next section. 



VI. NUMERICAL RESULTS 

In this section, we compare the heuristic and the an- 
alytical results of sections HH IIVI and [V] with the direct 
numerical solution of the full set of linearized equations: 



+ - = - - K ^ 2 *) . (109) 
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where £ = x/6cs, « = ^^cs and e = <5cs/£cs = 5 1/2 . 
These equations follow straightforwardly from Eqs. (|47l - 
I48|) after Fourier decomposing in the y-direction: Sip = 
exp(iKy), 50 = — i$(£) exp(iKy). The functions 
/(£)' ff(0; and w (0 are the normalized SP-like 

background equilibrium profiles: the reconnecting and 
reconnected magnetic field components and the in- 
flow and outflow velocity profiles, respectively (see sec- 
tion [TFf^] for the normalizations adopted for them). Ex- 
plicit expressions for these functions, parametrized by the 
position along the sheet yo = yo/Lcs, are derived in Ap- 
pendix [X] For the equilibrium adopted, the Alfvcn Mach 
point of the current sheet, defined by %(£ = 0, j/o.crit) = 
B y (£ = oo,y 0: crit), occurs at yo.crit ~ 0.61. Above this 
point, the upstream magnetic field is no longer able to 
stabilize the KH mode. 

Eqs. p09M110)l are solved in a domain of size 
—L x < £ < L x using a second-order-accurate predictor- 
corrector numerical scheme. The boundary conditions 
are *(-L Xi t) = $(L x ,t) = $(-L x ,t) = $(L x ,t) = 0. 
The size of the simulation domain L x depends on k, with 
lower values of k requiring larger domains [this is due 
to the behavior of the eigenfunction in the external re- 
gion, \I> ~ e _Ke ^; see Eq. (|60|) ]. Convergence tests were 
performed to ensure that both the domain size and reso- 
lution were appropriate. 



A. Plasmoid and KH instabilities 

We begin by focusing on the inviscid limit, Pm = 0. 
Plotted in Fig. [Ha) are solutions of: (i) the full analytical 
dispersion relation for the plasmoid instability, Eq. (|89p . 
for e = 10~ 6 (i.e., S = 10 12 ) and several different val- 
ues of yo (thin colored lines); (ii) the analytical KH dis- 
persion relation in the resistive limit, Eq. (|92[) (orange, 
long-dash-dot line), and in the ideal limit, Eq. (|107|) (red 
dotted line), for y = 0.8 (i.e., beyond the Alfven Mach 
point yo.crit = 0.61). Overplotted (symbols) are the nu- 
merical results. The theoretically predicted slopes for 
the plasmoid instability in the limits A< 1 and A — > 1~ 
[k 2 / 3 and K -2 / 5 , respectively; see Eqs. (|9~T1I9"2"1) ] are shown 
by the dotted black lines. The vertical dotted line is at 
ne = 1; all the analytical dispersion relations plotted arc 
only valid for values of k significantly to the left of this 
line. The vertical dash-dot-dot line shows the position of 
Ktr, the value of n at which a transition from the resistive 
to the ideal KH mode occurs, as given by Eq. (|108[) . 

For yo < 0.4, the agreement between the numerical 
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FIG. 2. (a) Growth rate (normalized to the global Alfven time ta = Lcs/Va) as a function of the wave number k — fcLcs for 
e = 10 -6 (i.e., S = 10 12 ) and at positions along the sheet, yo- Thin colored lines show the solution of the analytical dispersion 
relation for the plasmoid instability, Eq. (|89[) . The orange long-dash-dot line is the solution of the analytical resisitive KH 
dispersion relation, Eq. (|92p . and the red dotted line is the solution of the analytical ideal KH dispersion relation, Eq. (|107|l . 
both evaluated for yo = 0.8. Symbols are the results of direct numerical integration of the linearized equations. The dotted 
black lines identify the analytically predicted slopes for the plasmoid instability (k 2//3 for k < e -3 / 4 and k _2//5 for k 3> e -3 / 4 ). 
The vertical dotted line is at = 1; the vertical dash-dot-dot line identifies Kti (Eq. |TDSJ)). The analytical dispersion relations 
are numerically solved assuming that the global equilibrium can be approximately described by the analytical solution in the 
vicinity of the origin derived in Ref. [9^. For this equilibrium, yo.crit = 0.61. The numerical results use the analytical SP 
equilibrium calculated in Appendix El (for which yo,crit is the same), (b) Same as (a) for yo = 0.8 and e = 10 -8 (i.e., S = 10 16 ). 



solution and the analytical plasmoid dispersion relation 
is very good up to wave numbers approaching Ke ~ 1; 
this is not surprising since we used the ordering Kf < 1 
in our calculation. Importantly, this ordering is indeed 
respected by k = K m ax, Eq. (|93|) . which is accurately 
described by the analytical solution at these values of 
yo- Note that the dependence of 7 on yo is very weak, 
both as predicted by theory [Eq. with /q, w and /oo 
given by Eqs. (|A13|) . (|A"3]) and (|A14|) . respectively] and 
as determined by the numerical solution (see Fig. [4j). 

The behavior of the growth rate is distinctly differ- 
ent for y = 0.8. This value of yo is above the Alfven 
Mach point, yo.crit = 0.61, and it is, therefore, in the 
KH-unstable part of the current sheet. At low values of 
k, the analytical dispersion relation, Eq. (|92[) (labelled 
"KH res." in the figure) correctly captures the numeri- 
cal solution. The transition from the resistive KH to the 
ideal KH mode occurs at k « Ktr given by Eq. (11081) . 
The ideal KH dispersion relation derived in section [V] 
Eq. (|107|) (labelled "KH ideal" in the figure), applies for 
k > K tr , but fails to capture the fastest growing mode, 
which now occurs when ne ~ 1, as anticipated. 

The transition between the resistive and ideal modes 
is not completely clear in Fig. [2ja); for k > Ktr the nu- 
merical data points lie between the ideal and resistive 
KH lines (red, dotted, and orange, long-dash-dot, respec- 
tively) before rolling over and reaching 7 max at K = K max . 
This is because, even at such a large value of S, there 
is not enough scale separation between n tT ~ e -3 / 4 , and 
fcmax ~ e" 1 for the ideal solution to match its asymptotic 
form derived in section [Vj In order to illustrate clearly 



this transition, we plot in Fig. Hfb) the dispersion rela- 
tion obtained at yo = 0.8 for an even smaller (perhaps 
unrealistically so) e = 10~ 8 , i.e., S = 10 16 . In this figure, 
it is clear that the mode becomes ideal for k > Ktr, an d 
is correctly described there by the ideal KH dispersion 
relation, Eq. (fTUTl) . 

In Fig. [31 we plot the fastest growth rate and the cor- 
responding wave number as functions of the Lundquist 
number. Whereas the plasmoid instability scalings are 
obtained at y = 0.0 (i.e., 7 max - 5 1/4 , see Eqs. p i94]l ). 
we see that it is the KH scaling that is manifest at 
yo = 0.8 (i.e., 7 max ~ 5 1 / 2 , as derived in Eq. ([55])). 

Plots of 7 max and ft max as functions of yo at S = 10 12 
are shown in Fig. 2] In this figure, the dashed verti- 
cal line identifies the position of the Alfven Mach point, 
2/o,crit ~ 0.61. As expected based on the previous dis- 
cussion, the agreement between the analytical plasmoid 
dispersion relation, Eq. ([89]). and the numerical solution 
is excellent at values of yo < yo.crit- As yo — > yo,crit, the 
difference between the analytical and numerical solutions 
increases and explodes at yo = yo,crit- The transition to 
the KH mode happens then; for yo > yo,crit, our sim- 
plified analytical theory (section [Vj fails to produce a 
maximum of the growth rate and so cannot be compared 
to the numerical solution, which confirms that the rea- 
son for the failure of the asymptotic theory is that for 
7 = 7maxj K max£ is not small, but approaches values of 
order unity. 

Finally, in Fig. [5J we show the cigenfunctions "F^) 
(left) and $(£) (right) corresponding to the fastest grow- 
ing wavenumber, ft max , at S = 10 12 . These plots are con- 
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FIG. 3. Maximum growth rate (left) and the corresponding wave-number (right) as functions of the Lundquist number S, for 
j/o = 0.0 and yo = 0.8. The plasmoid instability [7 max ~ S 1 ^ 4 , K max ~ 5 3/ ' 8 ; see Eqs. (f9l I94|) ] observed at j7o = 0.0 is superseded 
by the KH instability [7 max ~ S 1/2 , K max ~ 5 1/2 ; see Eq. ([25])] at y = 0.8. 




FIG. 4. Maximum growth rate (left) and the corresponding wave-number (right) as functions of the position yo along the sheet, 
for e = 10~ 6 (i.e., 5" = 10 12 ). The analytic solution is given by Eqs. (|93H94[) . The vertical dotted line identifies the position of 
J/o,crit for the equilibrium parameters specified in Eq. (| A14[) . 



structed from runs at the values of yo plotted in Fig. @] 
We see that the eigcnfunctions undergo an abrupt change 
at yo ~ 0.6, where the Alfven Mach point of the system 
is located (identified by the dashed white line) . Close in- 
spection of the "J eigcnfunctions reveals that for j/o > 0.6, 
^(0) — > 0, as predicted in section [V] i.e., the most un- 
stable mode is ideal. Visible in the plot of the $ eigen- 
function beyond the Alfven Mach point is the formation 
of structure at each of the KH-unstablc shear layers, lo- 
cated at x/5qs ~ ±1. For a clearer observation of both 
these properties of the eigenfunctions, we plot in Fig. [6] 
one-dimensional cuts of Fig. [SJ taken at j/o = 0.4 (i.e., 
below the Alfven Mach point) and j/o = 0.8 (i.e., above 
the Alfven Mach point). As seen, for j/o = 0.4, ^ is 
finite at £ = 0, whereas it is zero at j/o = 0.8, in agree- 
ment with the analytical theory of section [V] and the 
prediction that above j/o,crit the fastest growing mode is 



non-reconnecting. Furthermore, the broadening of the \1/ 
cigenfunction around £ = at j/o = 0.8 suggests the pile- 
up of the magnetic field that we discussed in section [V] 



B. Effect of viscosity 

In order to address the effect of viscos- 
ity on the plasmoid instability, the term 
Pm ($"" - 2k 2 £ 2 $" + K 4 e 4 $) /k, where Pm = v/rf, 
is added to the right-hand side of Eq. (|110p . 
and two additional boundary conditions are used: 
&'(-L x ,t) = $"(L x ,t) = 0. We follow the generaliza- 
tion of the SP scalings to plasmas where Pm 1 derived 
in Ref. [13 : namely, we must scale the electric field at the 
origin as Eq/Ptti 1 / 4 , and the width of the current layer 
as 5qs — > SspPm 1 ^, where S S p/L C s = S^ 1 / 2 . There- 



15 




o 
i 




5. Eigcnfunctions (|\E'| 
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FIG. 

S = 10 12 ). These plots are constructed from runs at the values of yo shown in Fig. [4] In all runs, the size of the simulation 
domain is L x = 100, and the resolution is Ax = 0.0125. Only a fraction of the simulation domain is shown. The horizontal 
dashed white line shows the location of the Alfven Mach point, yo = yo,crit- 
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FIG. 6. Real and imaginary parts of the eigenfunctions for k 
Only a fraction of the simulation domain is shown. 



at e = 10~ 6 (i.e., S = 10 12 ) and y = 0.4 and yo = 0.8. 



fore, we rescale the parameter e = Scs/Lcs according to 
e e spP m1 ^ 4 where e$p = Ssp/Lcs = S^ 1 ^ 2 . 



Plotted in Fig. [7] are the maximum growth rate and 
the corresponding wave-number as a function of Pm, for 
S = 10 12 and yo = 0. Both K max and 7 max are seen to de- 
crease with increasing Prandtl number; a good fit to the 
data is given by 7 roax oc Pm -5 / 8 and ft max oc Pm~ 3 / 16 . 
The scaling of 7max and fv max with the Lundquist num- 
ber at Pm = 30 is shown in Fig. [5] We see that the 
S dependence of the maximum growth rate and of the 
corresponding wavenumber remains unchanged at large 
Pm, i.e., 7 m ax °c S 1 ^ and K max oc S 3 ^ s . These results 
agree exactly with the power laws derived in section [TTl 
Eqs. (HHHS1). 



VII. SUMMARY 



In this paper, a two-dimensional linear theory of the in- 
stability of large-aspect-ratio, Sweet-Parker-like current 
sheets is presented. This work is a direct generalization 
of our previous results [7(| (Paper I), where the simple 
equilibrium used was only a good model of a current sheet 
in the immediate vicinity of y = (y is the outflow di- 
rection). 

In the work presented here a general 2D SP-likc cur- 
rent sheet equilibrium is considered. As in Paper I, 
we conclude that large-aspect-ratio Sweet-Parker cur- 
rent sheets are violently unstable to high-wave-number 
tearing-like perturbations, and the same scalings of the 
growth rate with the Lundquist number S = LVa/tj are 
obtained here: 7 max TA ~ S 1/4 and k max Lcs ~ S 3 ^ 8 - 
see Eqs. (fMHiM)) . The plasmoid chain is formed inside a 
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FIG. 7. Maximum growth rate (left) and the corresponding wave-number (right) as a function of the Prandtl number Pm for 
e = 1(T 6 (i.e., S = 10 12 ) and y = 0. 
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FIG. 8. Maximum growth rate (left) and the corresponding wave-number (right) as a function of the Lundquist number S at 
Pm — 30 and yo — 0. 



boundary layer whose width scales as Si nner /8cs ~ S" 1 ^. 
These scalings have been confirmed via direct numerical 
simulation [43, [7?| ■ 

The more general approach employed in this paper has 
allowed us to calculate the growth rate of the plasmoid 
instability as a function of the position along the current 
sheet, yo- The dependence of 7 max and /c max on yo is a 
nontrivial function of the particular equilibrium consid- 
ered and, in the absence of a known analytical solution 
to the SP problem, cannot be evaluated explicitly. How- 
ever, for yo/Lcs *C 1 wc make use of the semi-analytical 
results of Uzdensky and Kulsrud [92[ and present an ex- 
act solution — Eqs. ([97l 158)1 . The most unstable wave- 
number and corresponding growth rate are then found 
to increase with distance from the center. Under general 
conditions (Syrovatskii-like upstream magnetic-field pro- 
file and outflow profile increasing monotonically along 
the layer), we show that the same result holds true at 
arbitrary yo /Lqs ~ 1 ■ This finding is somewhat counter- 



intuitive: a priori, one could expect that the increasing 
strength of the reconnected field along the sheet, as well 
as the shear in the ouflow (in the y direction), would 
provide a stabilising effect. Our calculation shows, how- 
ever, that both are irrelevant to the instability. An in- 
tuitive understanding of why that should be so can be 
gained by comparing the strength of the upstream and 
the downstream magnetic fields at the boundary of the 
inner (plasmoid) layer, £ = dinner: 



El 
B, 



Z~<5inner 



<5cs 



■_gl/2 _ £3/8 > x 



(111) 



i.e., even at the scale of the inner layer the reconnected 
field B x is completely overwhelmed by the reconnect- 
ing field By. The gradient of the background outflow 
in the y direction, whose length scale is ~ ^cSi is also 
unimportant because fcmax^cs 3> 1 everywhere in the 
sheet. At the periphery of the sheet, for yo > j/o,crit, 
where 2/o,crit/ics is equilibrium-dependent but otherwise 
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0(1), the current sheet becomes unstable to the Kelvin- 
Helmholtz (KH) instability driven by the velocity shear 
between the Alfvenic reconnection outflow and the sta- 
tionary upstream plasma. This occurs because the mag- 
nitude of the upstream magnetic field is a decreasing 
function of the outflow coordinate y and eventually be- 
comes smaller than the outflow speed (which is an in- 
creasing function ofy/Lcs)- At. and beyond, the Alfvcn 
Mach point, where this happens, the magnetic field can 
no longer stabilize the current sheet against the KH in- 
stability. 

We find that the KH instability of the sheet can ei- 
ther be resistive (i.e., induce reconnection at x = 0), or 
ideal (no reconnection), with lower values of /cLcs corre- 
sponding to the former, and larger values to the latter. 
The fastest growing KH mode, A: max Lcs ~ S 1 ^ 2 (i.e., 
^max^cs ^ 1), is an ideal, non-reconnecting mode. This 
is because reconnection cannot occur at the fast rates 
required by the fastest growing KH mode. A useful anal- 
ogy can be made with the Taylor (forced reconnection) 
problem (93j : since there are two shear layers, one on 
each side of the current sheet, the KH instability of the 
sheet is conceptually similar to a situation where pertur- 
bations at distant walls attempt to drive reconnection at 
a rational surface. In the Taylor problem, it is also found 
that the perturbations at the walls do not drive recon- 
nection in the initial stage. However, the same analogy 
suggests that as the ideal KH mode evolves into the non- 
linear regime, it will cause the upstream magnetic field to 
pile-up in the current layer, eventually leading to its re- 
connection. This KH-driven reconnection that occurs at 
Va > J/o.crit will give rise to a plasmoid chain, just as the 
"pure" plasmoid instability that is found at yo < j/o.crit- 
Therefore, in practice, it may be difficult to distinguish 
between the two situations. 

It is also worth noting that the basic KH instability 
mechanism that we have described here is completely 
general, i.e., it should apply to any reconnecting current 
sheet, not just to those that can be described by the re- 
duced MHD framework that we have adopted here: the 
only ingredient it requires is the existence of an Alfvenic 
Mach point somewhere along the layer. This should be 
a generic feature of most reconnecting current sheets. 
Whether the layer is collisional or collisionless may af- 
fect the dynamics of the KH instability, but its existence 
is not dependent on the plasma collisionality. In this 
respect, our findings may be related to recent observa- 
tions of the KH instability in collisionless simulations of 
guide-field reconnection [97j |. 

Finally, the effect of viscosity on the plasmoid instabil- 
ity has been addressed via numerical integration of the 
linearised set of equations. Our results are that in the 
limit Pm = vjr\ 1, the fastest growth rate and wave- 
number of the plasmoid instability scale as: 

7ma X ~ S^Pm-V* ~ L 1 ^ Vj/y/ 8 ^ 5 / 8 , (112) 
Kmax ~ S 3 / 8 Pm" 3 / 16 ~ L 3 J* vfV 3716 ^ 3716 - (H3) 



We have not performed a rigorous analytical calculation 
of the plasmoid instability in this limit, but we have 
been able to recover these scalings in a non-rigorous way 
from known results on the visco-tearing and resistive- 
kink modes [H|, via the rescaling of the background 
magnetic shear length a — > Sqs ~ LcsS^^Pm 1 /^ [3(| . 
Although these scalings are only expected to apply for 
S 3> Scrit , P m 3> 1 , where S CI u is the critical value of the 
Lundquist number for the current sheet to be plasmoid- 
unstable, they lead to the prediction that 

S clit ~ 10 4 Pm 1/2 , Pm > 1. (114) 

This result, as well as those of Eqs. (|112II113|) are con- 
crete predictions that can be tested in direct numerical 
simulations of MHD reconnection in the large magnetic 
Prandtl number regime. 
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Appendix A: Equilibrium Considerations 

Exact two-dimensional solutions of Eqs. (f2~Tl - |2"8"]) de- 
scribing a Sweet-Parker-like reconnecting current sheet 
are not known. In principle, these can be obtained 
by substituting the expressions for ip and <j) given by 
Eqs. (|5TMT|) into Eqs. (|2TH2"8"]l and equating equal pow- 
ers of (y - y )/L C s = y-yo- To lowest order in y- y , 
we obtain the following equations: 

«(£)/(£) - fovitMO = /'(£) - Eo, (Al) 

u(0v"(0 «(£)«"(£) = .9(0/" (0 - f(0g"(0 + 

Pm «'"(£), (A2) 

where we have used the normalizations of Eq. (|33p . ne- 
glected terms of order e 2 , and defined the normalized 
electric field E = L C sE /{B 2 S C s)- 
Evaluated at £ = 0, Eq. (fATj) yields: 

fo = E o ~ yo v o9o, (A3) 
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whereas for ^> 1 wc obtain from the same equation 



foo 



(A4) 



These expressions are exact; however, we see that 
Eqs. (|A1HA2[) are not a closed set, since there are only 
two equations and four unknowns: /(£), (the nor- 

malized reconnecting and reconnected magnetic field pro- 
files, respectively), and it(£), (the normalized inflow 
and outflow velocity profiles, respectively). This closure 
problem is introduced by the expansion in (y — yo) (recall 
the discussion of section Hn)) . In order to obtain a SP-likc 
equilibrium, which we require for our numerical solution, 
one has to close Eqs. (|A1HA2[) , e.g., by guessing two of 
the four unknown functions, and solving those equations 
for the other two. Any model of the equilibrium that can 
be found in this way is necessarily non- unique (i.e., de- 
pendent on the guesses required for the closure); however, 
we will see in what follows that a qualitatively satisfac- 
tory model of a SP current sheet can be obtained by this 
procedure. 

Let us introduce an auxiliary function, s(£), defined by 
the following equation: 



9® = -T<£) -*(£)■ 

Joo 

Then, from Eq. flAT) we obtain 
v(0 =^ 



(A5) 



± 



7£,* 2 (0 , /o„«(0/(0-/'(0 + ^> 



vl 



(A6) 



Eq. (|A2|) can also be easily solved in the limit Pm = 
(viscous effects in the equilibrium that we are about to 
derive can be modelled by a rescaling of the current sheet 
thickness, the outflow speed and the reconnection elec- 
tric field according to the SP relationships in the viscous 
regime derived in [36|). Using Eq. (|A"5j) . Eq. (|A"2|) be- 
comes 



no 

f(0 ' 



(A7) 



to be solved subject to the boundary conditions s(0) = 
UooVo/foo - go, where v = v(0) and g Q = g(0), and 
s'(0) = (we demand that both v(£) and g(£) are even 
functions). 

The general solution to this equation is 



*(0 = Oxf(0 + c 2 f($ 



(A8) 



[The lower limit of integration on the last term on the 
right-hand side of this expression need not be specified 
as C\ can be redefined to absorb the difference between 



different lower limits; note however that we take the lower 
limit to be finite, i.e., neither nor oo.] 

At this stage, the equilibrium problem is solved if we 
provide functional forms for the reconnecting magnetic 
field, /(C), and for the inflow velocity profile, u(£). The 
simplest choice for /(£) is the "Harris sheet" j98j : 



( -r-e 



A qualitatively plausible choice for u(£) is 

/(0 



u (0 = -Uc 



fo 



(A9) 



(A10) 



Substituting Eq. (|A9|) into Eq. (|A8|) and evaluating the 
integral explicitly, we obtain: 



= ( 90 ~ -^v 

JOO 



Joo \ Joo 



(All) 



where the constants of integration Ci, C2 have been 
chosen to satisfy the boundary conditions we specified 
for s(£). Substituting Eqs. jMEB into Eqs. (|A5f 
IA6|) yields explicit expressions for the two remaining un- 
knowns, the reconnected magnetic field and the out- 
flow velocity profile, v(£). Although it is not particularly 
enlightening to write down these expressions in explicit 
form, it is useful to evaluate <?(£) for £ 3> 1. It is 



g(0\ 



«>1 



± 



fc 



-v - .go 



A 

foo 



£ = ±gU- (A12) 



This expression is used in section PV Al to estimate the 
magnitude of g(£) for £ ^> 1. 

The last step in obtaining an analytical SP-likc equi- 
librium solution consists of determining Eq, vq, go and 
/oo, all of which can in principle be functions of j/o- A 
reasonable choice for is a Syrovatskii-like profile [86[ : 



foo = y 



1 



As for Eq, Vo, ffo, their values at yn = 
culated semi-analytically in Ref. (92f [9 

E (go = 0) = 1.075; 
go(yo = 0) = 0.642; 
v (y = 0) = 1.286. 



(A13) 



have been cal- 



(A14) 



The simplest choice is to assume that these values are 
constant along the sheet [note, though, that a linearly 
increasing dependence of the outflow and reconnected 
field profiles is already included in the normalizations, 
Eq. p}]. 

Examples of the equilibrium profiles obtained in this 
fashion are shown in Fig. [9l for yo = 0.4 (left) and 
y = 0.8 (right) (in Eq. (|A"6]) . the solution with the '+' 
sign is chosen). We see that these profiles retain all the 
qualitative features expected of a true SP equilibrium. 
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FIG. 9. Analytic SP-like equilibrium profiles (Eqs. (|A"5| . (g9j|, (fATO)) ') evaluated for t/ = 0.4 (left), and y = 0.8 (right). 

These equilibria are obtained by chosing the functional form of the upstream magnetic field, /(£)i an d imposing that the inflow 
profile be such that u(£) = — Uoo/(£)//oo. The lowest order (in (y — yo)/Lcs) Ohm's law and momentum equation can then be 
solved for the two remaining unknowns, namely the reconnected magnetic field, g(£), and the outflow, v(£). 



Note that for these parameters, the Alfven Mach point 
of the system occurs at j/o.crit = 0.61. 

The solution found here can be viewed as a generaliza- 
tion to the entire current sheet of the equilibrium derived 
by Biskamp [88[, which is only applicable for jjq = 0. As 
mentioned above, the equilibrium profiles obtained by 
this procedure, though exact, are not unique, since they 
depend on the guesses for /(£), u(£); another ansatz can, 
in principle, yield a different, but equally plausible, equi- 



librium. For the purposes of this paper, however, we do 
not believe this to be a serious constraint since we ex- 
pect both the plasmoid and the KH instabilities to be 
largely independent of the fine details of the background 
profiles; this certainly seems to be true for the plasmoid 
instability, as is suggested by the agreement between the 
theoretical predictions of Paper I using a very simplified 
equilibrium and subsequent numerical studies [43l l47j . 
The profiles we have derived are a convenient model for 
solving the linear problem, as we do in section IVT1 
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